Nature of the vortex-glass order in strongly type-II superconductors 
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The stability and the critical properties of the three-dimensional vortex-glass order in random 
type-II superconductors with point disorder is investigated in the unscreened limit based on a lattice 
XY model with a uniform field. By performing equilibrium Monte Carlo simulations for the system 
with periodic boundary conditions, the existence of a stable vortex-glass order is established in the 
unscreened limit. Estimated critical exponents are compared with those of the gauge-glass model. 
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In spite of extensive studies for a decade, the question 
of nature of the thermodynamic phase diagram of ran- 
dom high- Tc superconductors has remained unsettled. In 
zero field, the possibility of an exotic thermodynamic 
phase with broken time-reversal symmetry, called the 
chiral-glass phase, has been discussed [1]. Even more 
attention has been paid to the in-field properties. For 
superconductors with point disorder, possible existence 
of a thermodynamic phase called the vortex-glass (VG) 
phase, where the vortex is pinned on long length scale 
by randomly distributed point-pinning centers, was pro- 
posed [2]. In such a VG state, the phase of the con- 
densate wavefunction is frozen in time but randomly in 
space, with a vanishing linear resistivity p^. It is a truly 
superconducting state separated from the vortex-liquid 
phase with a nonzero pL via a continuous VG transition. 

Since cuprate high- Tc superconductors are extremely 
type-II superconductors where the London penetration 
depth A is much longer than the coherence length, it is 
important to clarify first whether the proposed VG state 
really exists in the type-II, unscreened limit A — > oo. 
Indeed, stability of the hypothetical VG state has been 
studied quite extensively by numerical model simulations 
[3-9] . Many have been based on a highly simplified model 
called the the gauge-glass model. Previous simulations 
on the three-dimensional (3D) gauge-glass model gave 
mutually consistent results that a continuous VG tran- 
sition occurred at a finite temperature characterized by 
the critical exponents, ~ 1.3, i] ~ —0.5, z ~ 4 — 5 [3-6]. 

The gauge-glass model has some drawbacks [3] . It is a 
spatially isotropic model without a net field threading the 
system, in contrast to the reality. Furthermore, source of 
quenched randomness is artificial. The gauge-glass model 
is a random fiux model where the quenched randomness 
occurs in the phase factor associated with the flux. In 
reality, the quenched component of the flux is uniform, 
nothing but the external field, and the quenched ran- 
domness occurs in the superconducting coupling or the 
pinning energy. It remains unclear whether these simpli- 
fications underlying the gauge-glass model really unaffect 
the basic physics of the VG ordering in 3D. 

Recently, several simulations were performed beyond 
the gauge-glass model [7-9]. The present author studied 
the type of the lattice XY model where the above lim- 



itations of the gauge-glass model were cured [7] . While 
the VG state was found to be stable, the estimated criti- 
cal exponents, particularly v ~ 2.2, differed significantly 
from those of the gauge-glass model, posing a possibility 
that the gauge-glass model lied in a different universality 
class. However, due to the effect of employed free bound- 
ary conditions, the estimated critical exponents might 
possibly be subject to large surface effect. Vestergren et 
al studied a random pinning model which took care of 
the above limitations of the gauge-glass model in a dif- 
ferent way, to obtain a finite-temperature VG transition 
characterized by the exponents, v ~ 0.7, z 2± 1.5 [8], 
which differed significantly from either those of Ref. [7] 
or from those of the gauge-glass model [4-6] . Olsson and 
Teitel claimed on the basis of their simulations on the 
lattice XY model with weak disorder that the VG order 
was not stable even in the unscreened limit [9]. Thus, 
once one tries to go beyond the gauge-glass model, the 
present theoretical situation seems quite confused. 

In the present paper, I study the lattice XY model of 
Ref. [7], but now with applying periodic boundary con- 
ditions, to overcome the finite-size effect originated from 
surface. The Hamiltonian considered is 
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where 9i is the phase of the condensate at the i-th site 
of a simple cubic lattice with N = sites, and the 
sum is taken over all nearest- neighbor pairs. Aij is a 
link variable associated with the vector potential due to 
uniform external magnetic field of intensity h applied in 
the z-direction. In the Landau gauge, it is given by 
K^={Af ,A\ ,Al) = {fd,hi^,Q), where I < ix < L denotes 
the x-coordinate of the site i. Quenched randomness 
occurs in the superconducting coupling which is as- 
sumed to be an independent random variable uniformly 
distributed between [0,2 J], J > being a typical cou- 
pling strength. I impose periodic boundary conditions in 
all directions in order to eliminate surface 'spins' which 
might contaminate the bulk critical behavior. The field 
intensity is chosen to be /i = 27r/4 (/ = 1/4). The lattice 
sizes are taken to be multiples of four, i.e., L = 8, 12, 16 
and 20. 

Simulation is performed based on the exchange MC 
method, where the systems at neighboring temperatures 
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arc occasionally exchanged [10]. Equilibration is checked 
by monitoring the stability of the results against at least 
three-times longer runs for a subset of samples. Sample 
average is taken over 980 (L = 8 and 12), 248 (L = 16) 
and 200 {L = 20) independent bond realizations. 

We run two independent sequences of systems (replica 

1 and 2) in parallel, and compute a complex overlap q 
between the local superconducting order parameters of 
the two replicas V"!^'^^ = <2xp(«6'j-^'^'). 
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where the summation is taken over all N = sites. In 
terms of the overlap q, the VG order parameter and the 
Binder ratio is calculated by 

,(^) = [<M^>], 9 = 2-11^, (3) 

where (• • •) represents the thermal average and [• • •] rep- 
resents the average over bond disorder. 

By inspecting the spatial pattern of the vortex snap- 
shot and by calculating the Fourier amplitude of the spa- 
tial distribution of vortices, it is checked that no periodic 

vortex-lattice order is formed in the system. 

The size and temperature dependence of the calculated 
Binder ratio is shown in Fig. 1(a). As can be seen from 
the figure, except for the smallest size L = 8, g for dif- 
ferent L > 12 cross at T/J = 0.82 ± 0.02, indicating that 
the VG transition occurs at a finite temperature. The 
data show a rather clear splay-out, in contrast to the 
near marginal merging behavior observed for the case of 
free boundary conditions [7]. In contrast to the case of 
the chiral-glass order in the XY spin glass [11], g does 
not exhibit a negative dip characteristic of one-step-like 
replica symmetry breaking. 

Further evidence of a finite-temperature transition is 
obtained from the mean-square currents [4] . The current 
in the jj. direction 1^ (/U = x, y, z) is defined by 

Iy. = jY^ Jijsm{ei-9j-Aij), (4) 
<u>ll/^ 

where the sum is taken over all nearest-neighbor bonds 
along the n direction. We calculate both the longitudinal 

(along the applied field) and the transverse (perpendicu- 
lar to the applied field) mean-square currents given by 

Il = [<h>% 4 = [<I.>' + <Iy>']. (5) 
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FIG. 1. Temperature and size dependence of (a) the bulk 
Binder ratio, and of (b) the transverse Binder ratio. 

As can be seen from Fig. 2 (a), the transverse mean- 
square current for L > 12 exhibits a crossing at T/J = 
0.81 ± 0.02, close to the crossing temperature of the 
Binder ratio. The longitudinal mean-square current also 
exhibits a crossing behavior, as is shown in Fig. 2(b). Ap- 
parently, finite-size effects are severer in /£ than in /|,, 
as the crossing for smaller sizes occurs at relative high 
temperature, which shifts down rapidly to lower temper- 
ature with increasing L, i.e., Tcross/ J > 1-1 {L = 8 and 
12), ~ 1.00 {L = 12 and 16), and ~ 0.89 (L = 16 and 
20) . Aside from this strong finite-size correction, the ob- 
served behavior seems consistent with the occurrence of 
a single bulk transition at T/J ~ 0.81 as suggested from 
the transverse mean-square current. 

In order to further confirm that both the transverse 
and longitudinal spatial components order simultane- 
ously, I also compute the transverse VG order parameter 
and the transverse Binder ratio [7]. These quantities are 
defined in terms of the layer-overlap gjj, for the fc-th xy- 
layer of the lattice, 

ii=^j:4''*4'\ (6) 

where i G k means the sum over all sites belonging to the 
fc-th layer, by 
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FIG. 2. Temperature and size dependence of (a) the trans- 
verse and (b) the longitudinal mean-square currents. 

The calculated is shown in Fig. 1(b). As can be 
seen from the figure, gr shows a behavior quite similar 
to g, exhibiting a clear crossing at T/J = 0.81 ± 0.02. 
This also indicates that the model exhibits a single bulk 
VG transition where both the transverse and longitudinal 
components order simultaneously. 

Next, we turn to the critical properties of the model. 
Standard finite-size scaling analysis performed for the 
bulk Binder ratio, with setting Tg/J = 0.82, yields 
u = 1.2±0.3, as shown in the inset of Fig. 3 where 
t = {T -Tg)/Tg. Then, from the scaling of the order 
parameter f/^\ the critical-point-decay exponent is esti- 
mated to be 77 = —0.5 ±0.1: see Fig. 3. The estimate 
of v is also corroborated by a finite-size scaling analysis 
of the transverse mean-square current which, with 
Tg/J = 0.81, yields = 1.0 ± 0.2. The corresponding 
finite-size scaling plot is given in Fig. 4. (Unfortunately, 
meaningful scaling analysis is not feasible for /£ due to 
the rapid shift of the crossing points with L.) In order 
to examine the possibility of anisotropic scaling, finite- 
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size scaling of gx and q>p (see the inset of Fig. 4) is also 
performed. 1 get = 1.2 ± 0.3 and r] = —0.5 ± 0.1, 
which agree within the errors with v and r] determined 
from g and g^^-*. No sign of anisotropic scaling is thus 
found. Combining all these estimates, I finally quote 
1.1 ±0.2, 77 = -0.5 ±0.1. 
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FIG. 3. Finite-size scaling plot of the VG order parame- 
ter (main panel) and of the bulk Binder ratio (inset), with 
Tg/J = 0.82, = 1.2 and 1 + ?7 = 0.5. 
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FIG. 4. Finite-size scaling plot of the transverse 
mean-square current (main panel) and of the transverse VG 
order parameter (inset), with Tg/J = 0.81. Best values of 
exponents are taken to he v = 1.0 for /|., while u = 1.2 and 
1 + 77 = 0.5 for q!^\ 



On comparing the exponent values obtained here with 
those obtained for the system with free boundary condi- 
tions [7], one sees that they differ considerably. In the 
range of lattice sizes studied here, the application of ei- 
ther periodic or free boundary significantly influences the 
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estimates of critical exponents. One naturaly expects 
that periodic boundary conditions give better estimates 
for the bulk exponents, since no surface 'spin' exists there 
which might contaminate the bulk critical properties. 

If one compares the present estimates of exponents 
with those of the the gauge-glass model (with periodic 
boundary), the values of both v and 77 are compatible 
with each other, i.e., v = 0.2 and r] = —0.5 ± 0.1 

for the present model vs. v = 1.39 ± 0.20 and r\ — 
—0.47 ± 0.07 for the gauge-glass model [6]. This coin- 
cidence might suggest that the present model belongs to 
the same universality class as the gauge-glass model. It 
thus appears that the differences in the form of quenched 
randomness and in the spatial anisotropy due to external 
fields are irrelevant to the critical properties of the VG 
transition. 

Finally, I wish to refer to the possible effect of screening 
in real VG ordering. In real type-II superconductors in- 
cluding high- Tc superconductors, the penetration depth 
is large, but of cource is not infinite. As several calcu- 
lations on the gauge-glass model and other random XY 
model have constantly suggested, finite screening effect 
might destabilize the stable VG state [7,12,13]. Since 
high-Tc superconductors are strongly type-II supercon- 
ductors, this screening- induced rounding of a sharp VG 
transition is a weak effect, visible only at very close to 
Tg of the unscreened system. With use of the present 
estimate of z^, a rough estimate of such a rounding (or 
crossover) temperature Tx may be obtained. Screen- 
ing effect would be visible when the coherence length 
C ~ (Co is the zero-temperature coherence length 

and t = {T — Tg)/Tg) grows comparable to the zero- 
temperature penetration depth Aq. Since the ratio Aq/^o 
is of order 10^ in high- Tc superconductors, the crossover 
temperature is of order ~ 10~^. It means that one 
has to approach Tg as close as t ~ 10~^ in order to see 
the screening-induced roimding. In other words, in the 
temperature range outside , the critical behavior of the 
unscreened system is expected to be observed experimen- 
tally [14]. 

In summary, the VG ordering of strongly type-II su- 
perconductors with point disorder is investigated based 
on a lattice XY model with a uniform field. The occur- 
rence of a finite-temperature VG transition is established 
in the unscreened limit. The estimated critical exponents 
V and T] are close to the corresponding gauge-glass val- 
ues, suggesting that the present model belong to the same 
universality class as the gauge-glass model. 

The author is thankful to Dr.S. Teitel for useful discus- 



sion. He is also indebted to Dr. J. Lidmar for his comment 
on the estimate of critical exponents. The numerical cal- 
culation was performed on the Hitachi SR8000 at the 
supercomputer center, ISSP, University of Tokyo. 

Note added 

In the course of preparing this manuscript, the author 
learned that Olsson also made a MC simulation of the 3D 
VG order in the unscreened limit based on a lattice XY 
model, with a diffierent choice of the random- coupling 
distribution [cond-mat/0301624]. Finite-temperature 
VG transition with an exponent j/ = 1.50 ±0.12 was ob- 
served, roughly being consistent with our present result. 
Critical exponents similar to our present values were aslo 
reported very recently by Lidmar based on a different 
type of VG model [cond-mat/0302577]. 
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